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. In this paper we investigate wormhole and spherically symmetric solutions in 4D gravity plus a 

matter source consisting of a ghost scalar field with a sine-Gordon potential. For the wormhole 
solutions we also include the possibility of electric and/or magnetic charges. For both types of 
solutions we perform a linear stability analysis and show that the wormhole solutions are stable 
and that when one turns on the electric and/or magnetic field the solution remains stable. The 
linear stability analysis of the spherically symmetric solutions indicates that they can be stable or 
unstable depending on one of the parameters of the system. This result for the spherically symmetric 
solution is nontrivial since a previous investigation of 4D gravity plus a ghost scalar field with a A^* 
interaction found only unstable spherically symmetric solutions. Both the wormhole and spherically 
I , symmetric solutions presented here asymptotically go to anti-de-Sitter space-time. 

I. INTRODUCTION 

(N 

> . 

' The discovery of the accelerated expansion of the present Universe, has led to a search for a possible mechanisms 
^sj , of this expansion. From the Einstein equations one finds that such an expansion is only possible if one violates 
one of the energy conditions. In hydrodynamical language, this means that the parameter of the equation of state 
T— I uu — - = — pressure should be less than —1/3 (i.e. violation of the strong energy condition), or even less than — 1 

• ^ p energy density ' ^ ^ ' ' 

\^ • (i.e. violation of the weak energy or null energy condition). In this last case where w < —1, one finds exponentially 

] fast expansion of the Universe. One of the simplest ways of obtaining such an acceleration is via models with ghost 

■ scalar fields. In these models, one assumes that the sign in front of the kinetic term of the scalar field is the opposite 
I ■ from that of the standard scalar field. It is well-known that ordinary scalar fields only lead to an equation of state 

^ with w > — 1. Ghost, scalar fields can give an equation of state with w < — 1 in which case they are called phantom 
fields [44] . Astronomical observations [ll, Q , indicate that violation of the weak energy condition is possible. Despite 
the strong theoretical misgivings about ghost fields these experimental observations have motivated the theoretical 

■ consideration of such ghost scalar fields. Depending on the choice of the potential energy of these ghost scalar fields 
one can obtain models of the accelerated expansion of the Universe that are in good agreement with the observational 
data 0,13. 

Ghost fields have also been studied in the context of compact astrophysical objects, such as stars and wormholes 
composed partly or completely of ghost fields. In particular, ghost fields find a natural application to wormholes 
since these objects are know to generically require violation of the weak (or null) energy condition [5|. In the papers 
[3 ^ ghost scalar field with a Mexican hat potential was investigated and it was found that this system only had 
regular, stable solutions for topologically non-trivial (wormhole-like) geometry. The works [8l-[l3| further investigated 
traversable Lorentzian wormholes refining the conditions on the type of matter or field that would lead to such space- 
times. In the papers [3, [isj wormhole solutions induced by quantum effects in GUTS were proposed. A general 
overview of the subject of Lorentzian wormholes and violations of the various energy conditions can be found in the 
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book by Visser [T^ . 

There have also been papers devoted to non-wormhole hke solutions which have a trivial topology. The paper 
[13] gives a general classification of nonsingular, static spherically symmetric solutions to the Einstein equations for 
scalar fields with arbitrary potentials and negative kinetic energy. In Ts^ it was found that there existed spherically 
symmetric solutions to the system of 4D gravity plus a ghost field with a Mexican hat potential which were similar 
to the Bartnik-McKinnion solutions [l^ but with the ghost field replacing the Yang-Mills field. A linear stability 
analysis of the solutions found in [Tsj indicated that these solutions were not stable, although depending on the form 
of the effective potential they might persist for long times - i.e. they could be meta-stable. 

Having a negative kinetic energy term leads to obvious quantum instability, but there are possible solutions of this 
problem of quantum instability [20,[2l|. However, the question of the classical, linear stability of a solution with ghost 
scalar fields requires a case by case study, for example: (i) In Q it is shown that the solutions with the Mexican hat 
potential are stable against radial perturbations, (ii) In [2^ the classical, stability conditions for a charged wormhole 
created by a massless, ghost scalar field were found, (iii) In [2^ [23] it is shown that wormhole solutions within the 
general framework of scalar-tensor theories with massless, non-minimally coupled scalar fields are unstable under 
spherically symmetric perturbations, (iv) In [2^, [2^ the generalization of this problem was carried out for the case 
with electric or magnetic charges, (v) In [13] a linear stability analysis for a wormhole created by a massless ghost 
scalar field was carried out with the result that these wormholes had two asymptotically spherically symmetric ends 
with cither both ends having vanishing ADM masses or with the ends having ADM masses of opposite signs. It was 
found that such wormholes are stable if the ADM masses are sufficiently small, (vi) In the paper 27] some estimates 
of the stability of a rotating, scalar field wormhole were given, (vii) In [2^ [2^ it was shown that wormholes with a 
massless ghost scalar field are unstable with respect to linear and nonlinear perturbations, (viii) A finally example of 
such classical, stability analysis can be found in [s^] where it is shown that adding a charge to the model of (vii) did 
not lead to a stable wormhole solution. 

These examples show that the question of classical stability or instability of a spherically symmetric or wormhole 
configuration is strongly mo del- dependent and that one needs to carry out a stability analysis in each specific case. 
In this paper, we carry out a linear stability analysis for a ghost scalar field model with the a sine-Gordon potential. 
We also show why this system does not admit a topological or global stability analysis. In a previous paper [3l| this 
system was studied with the interesting result that one could use the solutions to generate a dynamical model for the 
creation/annihilation of wormholes - by allowing the throat radius to vary smoothly from some finite value down to 
zero one could interpret this process as the creation a wormhole with subsequent annihilation into two, disconnected 
spherically symmetric solutions. In this paper, we will generalize the solutions of [Slj to the include electric and/or 
magnetic charges, and carry out a linear stability analysis for all these solutions. 

For the solutions found in [3l| it was shown that the static spherically symmetric solutions with a ghost field plus 
sine-Gordon potential tended asymptotically to anti-de-Sitter (AdS) spacetime. AdS space-times have generated much 
interest due to the discovery of the AdS/CFT correspondence [33H35i] which implies that a gravitational theory in an 
AdS space-time is dual to a conformal field theory (CFT) on the boundary of the space-time. The question of the 
stability of such space-times has already received considerable attention [3a - [38| . In [sg , in the context of Einstein's 
theory plus a negative cosmological constant, the stability of asymptotically AdS metrics against linear perturbations 
was demonstrated. In contrast reference [37| pointed out that in some cases instabilities may arise for asymptotically 
AdS spaces. This implies the necessity of performing a stability analysis for solutions with asymptotically AdS 
space-times on a case by case basis. 

A second motivation to study the stability of the solutions presented below is related to the existence of the stable, 
static, sine-Gordon solitonic solutions for ordinary scalar field in (1-1-1) dimensions without gravity [39|. The sine- 
Gordon solutions are stable because of topological reasons. In this paper we show that for some range of parameters 
adding gravity and changing the ordinary scalar field to a ghost scalar field leads to classically, linear stable solutions 
in (3-1-1) dimensions. For other ranges of the parameters the solutions are unstable. Thus including gravity and 
allowing ghost fields leads to linear stable solutions which is in contrast to ordinary scalar fields in (3-1-1) which do 
not have stable solutions. However, in contrast to the sine-Gordon soliton kink solution, the present solutions do not 
obtain their stability from topological arguments. 



II. FIELD EQUATIONS WITH ELECTRIC AND MAGNETIC CHARGES 

We consider a model of a gravitating ghost scalar field in the presence of the electromagnetic field Fik- The 
Lagrangian for this system is [3l| 

L = -ji^- - V{^) - , (1) 



3 



where R is the scalar curvature, G is Newton's gravitational constant, and V is the sine-Gordon potential with a 
reversed sign from the usual case 



V 



m 
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Here m is a mass of the field, and A is a coupling constant. The corresponding energy-momentum tensor is 



(2) 



(3) 



To perform a linear stability analysis for this system, we take a time dependent, spherically symmetric metric of 
the form 



A{r,t) 



(4) 



where the metric functions F{r,t) and A{r,t) depend both on the radial coordinate, r, and time, t. 

Let us consider only radial components of electric and magnetic fields Fqi — Er and F23 = —Hr- Then by choosing 
the ansatz for the magnetic field F23 = dgA^, where A3 = Qcos9 and Q is the magnetic charge, we have the following 
expression for the quadratic field combination 



-2Ae'^^E^ + 2 



(r2 -t-r2)2' 



The equation for Er can be found from Maxwell's equations 

[V^i^'1.,=0 Er = 



{r'^+r^)VA' 



(5) 
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where q is the electric charge. 

Using the Lagrangian ([T|) and equations (O and ([HI), the gravitational equations take the form 
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We have introduced new dimensionless variables (/) — {\/\/m)ip, x — mr,T = mt and j3 — SnGm^/X, q = ^/Xq,Q 
VXQ. Equations ^ and dH) are the (l) and [(Jj.) — (^)] components of the Einstein equations. 
The scalar field equation is 
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In the above equations, "prime" and "dot" denote differentiation with respect to x and r, respectively. 

In [3l| a set of static solutions was found for the system (I?])-® for various values of the parameter xq and with the 
electric and magnetic charges (g, Q) set to zero. It was shown that there exist regular solutions up to the value xq = 0. 
This was interpreted as the possible creation/annihilation process of a wormhole - as the throat radius a:o — >■ the 
wormhole split into two, separate, spherically symmetric space-times. This interpretation assumed that the static 
solutions could be considered as snapshots of a dynamical process where the wormhole annihilated into two separate 
space-times or run in reverse where two separate space-times merged to create a wormhole. These solutions and this 
creation/annihilation interpretation are similar to the 5D Kaluza-Klein dyonic wormholes studied in 32 | 

We now extend the results of j3l| by including radial electric and magnetic fields. These fields will give positive 
contributions to the total energy density of the system, leading to solutions differing from the ones discussed in . 
We have not found an analytical solution for the system ([T])-® for the entire range of x. Thus we solve this system 
of coupled, non- linear equations numerically. From the equations (O-® we obtain the initial conditions for starting 
our numerical solution as 



F(0) = const., ^(0) = 1-/3 
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(10) 
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FIG. 1: The background, static, wormhole solutions of the 
system (l?])-© for the scalar field </) at xo = 0.3, /? = 1 and 
g = (the solid line), q = 0.35 (the dashed line), Q = for all 
graphs. Asymptotically (f> ^ n. 
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FIG. 2: The background, static, wormhole solutions of the 
system ((Tl)-® for the dimensionless electric field Er — 

{^/\/m^)Er = qe^ I [(x^ + xg) V^] at = 0.3, /? = 1 and 

q = 0.35, q = 0.40, and q = 0.43, top-down, respectively, 
Q = for all graphs. Asymptotically E — >■ q/x'^. 



These initial conditions, where xq 0, give wormhole solutions. Later we will discuss the spherically symmetric, 
boson star sokitions for which xq = 0. From the above expression for ^4(0) one can see that the addition of a nonzero 
electric and/or magnetic charges q, Q gives a lower limit on xq for which A{0) is positive. We do not give the explicit 
expression for xq but it involves solving a quadratic equation for Xq. The physical meaning is that the addition of a 
large charge - and its associated positive energy density - makes the existence of wormholes and spherically symmetric 
solutions impossible. The general, qualitative behavior of the solutions at acceptable values of xq (i.e. values for which 
A{0) > 0) with electric and/or magnetic fields remains similar to the case when there are no fields (compare figures 
[T]|3]of the present paper with figures 1-3 in [3l|). From now on we take into account only electric charge by taking 
q non-zero and setting Q = 0. This is not a reduction in the generality of the system since, from equation (O, it is 
possible to introduce a new effective electric charge defined as g^yy = (f + Q^- As implied by electromagnetic duality 
it is only a matter of convention whether one talks of electric charge, magnetic charge or some combination. What 
matters is qef / . From the numerical solutions shown in figures [T][3] one can make the following statements about the 
asymptotic form of the solutions: (i) The scalar field approaches the minimum of the potential V{(p) from ^ (i.e. 
(j) = it) both from below (for a; < 0) and from above (for x > 0). (ii) The electric/magnetic field asymptotically takes 
a Coulomb form, I/x^. Note as the effective electric/magnetic charge increases a cusp develops a,t x — 0. This might 
be seen as the influence of the electric/magnetic charge to generate a singularity at a; = in the case of a space-time 
with trivial topology i.e. when — > and the wormhole space-time goes over to the spherically symmetric topology, 
(iii) The metric ansatz functions F{x),A{x) asymptotically go as e^^ — > ^x'^ and A — ^2;^ i.e. asymptotically the 
metric approaches Anti-de-Sitter space-time. In order to have the space-time be asymptotically AdS one needs to 
choose the constant value of F{0) = const, such that both e^^^^^^ and A{x) approach ^x'^ as x 00. In figure[31 for 
example, F{0) « —0.945 at q = 0.35 and F{0) « 0.145 at g = in order to meet this requirement. The choice of F{0) 
is connected with the scaling of the time coordinate at x — 0. Thus one must appropriately scale the time coordinate 
at a; = in order to have both metric functions, e^^^^^ and A{x), give asymptotic AdS space-time. 

III. LINEAR STABILITY ANALYSIS 

We now study the dynamical stability of the above static solutions under linear perturbations with a harmonic time 
dependence. This type of stability analysis is essentially that used in the works [13, |4lj to study boson stars. The 
present authors have used this type of stability analysis to study the issue of stability for the spherical solutions for 
the system of a ghost scalar field with a Mexican hat potential coupled to gravity 18]. For simplicity, we introduce a 
new metric function a{x) in place of A{x) defined as A{x) = e°^^^. We perturb the solutions of the system (l7])-(l9|) by 



FIG. 3: The background static, wormhole solutions of the system (0-((9]) for the metric functions F and a = \nA at xo — 0.3, 
l3 = 1 and q = (the sohd lines), q = 0.35 (the dashed lines), Q = for all graphs. Asymptotically e^^ — > and A — >■ ^x^. 

This means that the space-time for large x is Anti-de-Sitter. 



expanding the metric functions and scalar field function to first order as follows 

F{x,t) = Fo{x) + Fi{x) cosujT, 
a(a;, r) — aQ{x) + ai{x) cos lot, 

(I){x,t) = 00 (2^) + 01(2;) 



\/a;2 + xl 



The index indicates the static background solutions of equations ([21)- (O and r is the scaled time parameter defined 
in the previous section. The first-order perturbation equations which follow from l^i and ([8]) are 
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and the equation for (j>i is 
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with the potential given by 

Vo{x)=e' 

Introducing the new independent variable p 
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we can rewrite equation (jl3p in a Schrodinger-like form 

~^ + V[x{pM,^u;'<I>,, 



(11) 
(12) 

(13) 

(14) 

(15) 

(16) 
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FIG. 4: The potential from (|17p as a function of p which is defined in (|15|l . The sohd line corresponds to xo — 1, the 

dash-dotted, the long dashed and the short dashed lines correspond to xq — 0.5, xq = 0.1 and xq = 0.01, respectively. The 
curve for a;o = (the spherically symmetric solution) is also shown by the solid line. For the spherically symmetric solution 
the range of p is < p < oo. For all solutions the electric and magnetic charges are set to zero (g, Q — 0) and 13 = 1. 



where 

V[x{p)] = e'^<'+''oVoix). (17) 

The study of the stabihty of the static solutions displayed in figures [T][3] now reduces to a study of the Schrodinger-Uke 
equation given by and ^T7\ . If the energy-hke eigenvalues, from (fTB|) are > then the solution is stable; if 
there are eigenvalues such that cj^ < then the perturbations grow exponentially and the solution is unstable. 

Before moving on to the details of the linear stability analysis of the wormhole and spherically symmetric cases 
we discuss several other criteria for stability of solutions similar to those presented here. There are several other 



stability criteria - topological stability [39| and global stability (42| - which are more powerful and go beyond the 
linear stability analysis presented here. However the present solution does not meet the conditions that allow us to 
apply these broader stability criteria. 

The typical example of topological stability is given by the sine-Gordon kink in the absence of gravity The 
sine-Gordon kink is topologically stable since the value of its scalar field takes two different values at a; = —oo and 
X — +00. From figure [T] one can see that the wormhole solution does not have this feature; the scalar field of the 
wormhole solution approaches the same value at a; = — oo and x = +oo thus ruling out a topological argument its 
stability. 

A global stability analysis method for systems of gravity coupled to a complex scalar field was given in (4^ . In this 
work the global stability was determined by studying the binding energy of the system 

B = M - mN 

where M is that Tolman mass and N is the particle number. The Tolman mass is defined by 

M = j{2T^~Ti:)^\Sx . 

From ([3]) — dfj_(j)d'^(t) + 4^(0). Combining this with Tq , which also can be obtained from (|3]) or read off of ([7]), we 
find that (2Tq — T^) in the above equation tends asymptotically to a constant value equal to — 2T4o = 4, where Voo 
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FIG. 5: The dependence of the lowest eigenvalues cj^ on the throat size xo without electric charge {q = 0) and with q — 0.3, 
P = 1. For q — 0.3 at xo ~ 0.21, tj^ — >■ and A(0) — >■ simultaneously (see equation (|10p 'l. 



is the asymptotic value of the potential. The asymptotic value of ^/\g\ oc and so the integral for M above diverges 
- as one would expect for a space-time which is asymptotically AdS. 

Next, the particle number N is the integral over all space of the time component, of the Noether current density 
J'' = ^V\9\9^''[rducjj~cbd,r] 

N ^ j fd?x . 

Since our scalar field is not complex — and so = 0. Thus neither the wormhole nor spherically symmetric 
solutions presented in section II can be analyzed in terms of topological stability or global stability due to the 
asymptotic AdS nature of the space-time, the asymptotic behavior of the scalar field or the fact that the scalar field 
is real rather than complex. 

A. Wormhole case (a;o 7^ 0) 

We now study the stability of the wormhole solutions for various values of the wormhole throat size xq. Figure 
|4] shows four effective potential curves for a;o = l,xo = 0.5, = 0.1 and xq — 0.01 respectively and with no 
electromagnetic field. One can see that for xq = 1 and Xq = 0.5 the curves of the potential are completely positive so 
that oj^ > for any solutions of ([T6| and (fT7| i.e. the solutions are stable. In the cases when xo = 0.1 and xq = 0.01 
there are small regions where the potential (jl7p becomes negative. From figure |4] these negative regions occur around 
p « 1.4. Thus in principle the lowest eigenvalue could be negative. Further analytical details which support 
the positiveness of for the parameters used in this paper are given in appendix A, even in those cases when the 
potential has regions where it dips below zero. In the language of ordinary quantum mechanics these regions 

are not deep enough for a bound state with cj^ < to form. 

For comparison figure |4] also shows the behavior of the potential ([T7| for the spherically symmetric solution, i.e., 
when Xq = and /3 = 1. The full investigation of this case is given in the next subsection and appendix B. One can 
see that near p = behavior of the spherically symmetric case differs from the behavior of the wormhole solution 
as Xq — 0. This difference is connected with different behavior of the derivative 4>'{Q) at x = 0: For the wormhole 
solution, this derivative tends to infinity, according to (|10p ; for the spherically symmetric solution, the derivative goes 
to [3l|. Nevertheless, at /3 = 1, the spherically symmetric solution is stable, since the lowest eigenvalue, w^, is 
positive. 
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FIG. 6: The spherically symmetric case: The scalar field (j) at 
different values of the parameter /3: 1.5, 1.0, 0.75 and 0.5, from 
top to bottom. Asymptotically — >■ vr. 




FIG. 7: The spherically symmetric case: The metric functions 
F{x) (the dashed lines) and a — In A (the solid lines) at dif- 
ferent values of the parameter /3: 1.5, 1.0 and 0.5, from top to 
bottom. Asymptotically e^^ and A 



3 ^ 



B. Spherically symmetric case {xq = 0) 



In this subsection, we study the stability of the system of equations ([Zl)-® when the radius of the wormhole throat 
goes to zero xo — > 0, i.e., the spherically symmetric case. The static solutions for this case were found in [3l| as a 
limiting case of the wormhole solutions when — > 0. In this case, the only free parameter is (3. We specify the 
boundary conditions as in [sij 



ylo(O) = 1, 00 (0) = const., 0;,(O) = 0, i^o(O) = const. 



(18) 



With these initial conditions wc numerically found solutions for (t>{x), F{x), A(x) for different values of the parameter 
/3 - see figures [6l [7l The boundary conditions from ()18|) also imply we are focusing on the case where the charges 
are set to zero since in the presence of charges the solutions become singular at x — ^ 0. This can be seen from the 
boundary condition for A{x) implied by equation (O: 



Ao(0) 



1 



2 



Some further analytical details of the linear stability analysis for the spherically symmetric case are given in appendix 
B. The summary of the results from appendix B are that for the spherically symmetric (xq = 0) case, the solutions 
are stable for /3 > 3/4 and unstable for P < 3/4. 



IV. CONCLUSION 



In the paper we have studied dynamical stability of spherically symmetric and wormhole configurations for the 
system of 4D gravity coupled to a sine-Gordon, ghost scalar field. In some we included cases electric and/or magnetic 
charges. The stability of these solutions was tested with respect to linear, radial, harmonic perturbations. For the 
wormhole configurations we considered the case both with and without electric and/or magnetic charges. For the 
spherically symmetric case we investigated only the case when the electric and magnetic charges were zero, since 
adding charges to the spherically symmetric solution would result in a singularity at a; = 0. 

To determine whether a particular configuration was stable or not the equation for the perturbation of the scalar 
field, ([13]), was written in a Schrodinger-like form The "energy" eigenvalue of this equation was the square of 

the frequency, w^, of the harmonic perturbations. If > then the solution was stable; if < this implied 
exponential growth of the perturbation and that the solution was unstable. This method of testing for stability was 
used in [13, \^ to study the stability of boson stars, and was used by the present authors fl8\ to study the stability 
of phantom/ghost boson stars. Other, stronger stability criteria - such as topological stability fsoj or the global 
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stability of f42| - were not applicable to the present solution due to the asymptotic AdS nature of the space-time, the 
asymptotic nature of the scalar field, or the real character of the scalar field. 

For the wormhole solutions, using a combination of numerical calculations and analytical estimates, it was shown 
that stable solutions exist for a wide range of the parameters (e.g. the size of the wormhole throat xq and the free 
parameter /3) both with and without the charges. For the spherically symmetric solution (when the wormhole throat 
radius xq =0), it is shown that the stability of the system depends on the value of the parameter /3: when /? > 3/4 
the solutions were stable; when /3 < 3/4 the solutions were unstable. 

The stability of these solutions, especially the spherically symmetric solutions, is a non-trivial result. Previously 
the authors [1^ found spherically symmetric solutions to the system of 4D gravity plus a ghost scalar field having 
a At/)'' self interaction. These solutions were compared to the Bartnik-McKinnon solutions [l^ with the ghost field 
playing the role of the SU(2) Yang-Mills fields. However these ghost field supported solutions were not stable. Here 
we find stable solutions by considering a ghost field with a sine-Gordon potential. In (1 + 1) there exists the well 
known sine-Gordon soliton [39| but it is known from Derrick's theorem [l^l that soliton solutions coming from only 
scalar fields do not exist in dimensions higher than (1 + 1). What we have found here is a soliton- like solution in 
(3 -I- 1) dimensions. This was achieved by adding gravity and letting the scalar field be a ghost field. The only 
feature about the spherical solutions that is slightly undesirable is that fact that asymptotically the solution goes to 
AdS space-time rather than Minkowski space-time. However, in light of the AdS/CFT correspondence [s^ and the 
connection between AdS space-time and the holography ideas jssf , this feature of the asymptotic form of the solution 
may not be so negative. Finally, if the universe has some cosmological constant as observations indicate, having a 
space-time with asymptotic dS or AdS behavior is actually preferable to having asymptotic Minkowski behavior. 
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Appendix A: Details of Wormhole Case (xq ^ 0) 

In this appendix we give additional, analytical support for the positiveness of the potential ([T7| (and therefore of 
w^) by analyzing the asymptotic behavior of l/[a;(/9)] for x ^ and for x — > oo. First we consider the behavior of 
the potential (IT71) at a; 0, making use of the boundary conditions pUj) . To simplify the analysis we express the 
condition for the metric function A (i.e. the middle equation in (|10p ) in term of a — In A. This yields 

ao{0)=ln{l-/3Sxl), (19) 

where the index refers to a static solution, and we have introduced the notation 

^=-4-2, 7= 5 >0. 

It follows from (fT9|) that A{0) > 0. By taking the limit A{0) = one finds a lower bound on xq, which we denote as 
the critical value, Xcr- 



2 _ 2 ^ 1 

''O ^cr > ^ 



-l + y/l-H 8/3^7 



(20) 



Taking into account this critical value, the behavior of the potential pT)) near x = depends on whether the parameter 
5 is positive or negative. In the case when (5 < its value falls in the range 

-2 < (5 < . 

The lower bound, 5 = —2, corresponds to absence of any charges (7 = 0). There are two limiting cases: 
(i) Xq ^ 1, /? ~ 1: in this case 

ao(0) « 0, ^^o(O) « 0, <^o « 0o(O) + (j^'^x 
and, accordingly, the potential pT)) will then be 

^[0] « 1 + 4 » 1- 
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In other words the potential wih be positive and large near a; = 0. 

(ii) Xq ^ 1: in this case ao(0) « In (/3|(5|2;q) and the potential P7|) will be 

V[0]^\S\ (1-+I3 



Again the potential is positive near x = 0. 

We now turn to the asymptotic behavior of the potential ((T7)) as x — cxd. In this case, the metric functions, 
F{x), A{x) behave as (see the discussion at the end of section II and the caption for figure [3]): 

e^^° -> ^x^, ao In 

To investigate the x — > oo behavior of the scalar field </> we expand the field as a small perturbation, S(l)oo , around the 
asymptotic solution of tt 

= 7r + (50oo, (5(/>oo < 1- 

Plugging this form into equation ([9]) yields the following equation for Scfioo 

4 3 
X Ix'^ 

This equation has the solution 

where C is an integration constant. Taking into account these asymptotic forms for the metric functions and the 
scalar field, the behavior of the potential p7|) is 

V[x{p)U « 

In summary for the case 5 < Q wc have the following behavior of the potential from \n\ near x — Q and as x oo 
f (l + ^) > 1 at xo < 1, /3 « 1; 
1 \5\ {\+p)>0 at xo > 1; 



X ^ 0: V[x[p)] 



X — oo: T^[x(/9)] -> gX^, xo,/3 are arbitrary. 

For the case 5 > 0, there is an upper limit, 5cr, coming from (j20p . In particular, at 7 — >■ 00, 5cr — >■ 0, and at 7 — >■ 0, 
5cr ^ l/(/3^7) ^ 00. Carrying out an analysis similar to the 5 <Q case, we find the following asymptotic behavior 
for the potential from p7|) for < (5 < 5cr'- 

> at (5 — > 0, xo > Xcr and finite; 




X ^ 0: V[x{p)] 

X — )• CX): y[x(p)] |x^, Xo > Xcr- 

The above analysis shows that near x = and at x — ^ 00 the potential (fT7|) is positive. Although we have no similar 
demonstration for the positiveness of the potential at intermediate x, the numerical solutions of l^[x(p)] show that for 
some of the parameters we investigated the potential was positive over the entire range of x - see figures SI [S] This 
shows that at least for these parameters the solutions found in section II and in are stable. 

Appendix B: Details of Spherical Symmetric Case (xq — 0) 

As for the wormhole solutions, we now study the stability of the static solutions by looking at the behavior of the 
potential (|17p near x — and x — > 00. First we expand the functions (/iq, Aq, Fq in a Taylor series in the neighborhood 
of X = 

0o«<^o(O) + 02y, (21) 
Ao«Ao(0) + A2y, (22) 
Fo«-Fo(0) + F2y, (23) 
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where 02,^2,^2 are the corresponding values of the second derivatives at x = 0. Expressions for these second 
derivatives can be found by inserting into the static equations ©-dH])- This yields 

02 = ^ sin 0(0), ^2 -/3[l-cos(/.(0)], F2 = ^A2. 
Inserting the expansions from ((2T|) - (|23| into the potential pT)) . gives 

V[0] = e2-P^°(°) [/3 + (1 - /3) cos ^{0)] . (24) 



1 1 1 1 1 1 1 1 1 1 1 1 

■ V 

1 
/ 
t 


i-T 1 1 1 


-n 1 1 1 y 


1 , , , , 1 , , , , 1 


i 

*\ 

*\ 

\ 

\ 

\ 



FIG. 8: The potential V[x{p)] from Hill as a function of p, where p is defined in The the values = 1.5, 1.0,0.76,0.74 

correspond respectively to the short dashed line, the solid line, the long dashed line, and the dash-dotted line. The critical 
value corresponds to P — 0.75, see Eq. (|25p . 



We now turn to the x — cx) behavior of the potential p7)) . In this case, the metric functions behave as 



As for the wormhole solutions, in order to study the behavior of the scalar field 4> at infinity, we take a perturbation 
around the asymptotic solution of the static equation ([9]) of the form 



With this equation ^ takes the form 



If /5 ^ I this has the following solution 



< 1. 



where C is an integration constant. If /? < I the solution is 



sin(^i^|Hlna;) cos(^^2lnx 



1 2.3/2 2 X'^/^ 
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where Ci,C2 are integration constants. 

Using these asymptotic expressions for the metric functions and the scalar field, we can obtain the asymptotic 
behavior of the potential (|17p for the two cases: 

(i) For /3 > 2/3: 



where 



— 1.5 < n < since 



^ 2/3 + 2"^ ^ 



^ Cx-^-^ at /3 ^ |; 
C'a;-i/(2/3) at ^ 00. 



(ii) For ^ < 2/3: 



2/3 



where the notation [cos, sin] denotes a set of terms containing the functions 

and sin 



cos 



27? 



In:; 



In; 



Hence, one can see that, asymptotically, the terms proportional to x^" in the first case and l/x in the second case 
can be neglected. The final expression for the asymptotic form of the potential for both cases (/? < 2/3 and /3 > 2/3) 
is 



V[x{p)] 



(25) 



It is clear from this expression that there is critical value of /? = 3/4. For /3 < 3/4 the potential becomes negative and 
the static solutions of the system ©-(O are unstable. 

A graphic illustration of the behavior of the potential ([T7| at various values of /3 is presented in figure [3 and the 
behavior of the lowest eigenvalues uP' as the dependence on (5 is shown in figure [SI 
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